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We identify ‘minimal seeds’ for turbulence, i.e. initial conditions of the smallest possible 
total perturbation energy density Ec that trigger turbulence from the laminar state, in 
stratified plane Couette flow, the flow between two horizontal plates of separation 2H, 
moving with relative velocity 2AU, across which a constant density difference 2Ap from 
a reference density is maintained. To find minimal seeds, we use the ‘direct-adjoint- 
looping’ (DAL) method for finding nonlinear optimal perturbations that optimise the 
time averaged total dissipation of energy in the flow. These minimal seeds are located ad¬ 
jacent to the edge manifold, the manifold in state space that separates trajectories which 
transition to turbulence from those which eventually decay to the laminar state. The edge 
manifold is also the stable manifold of the system’s ‘edge state’. Therefore, the trajecto¬ 
ries from the minimal seed initial conditions spend a large amount of time in the vicinity 
of some states: the edge state; another state contained within the edge manifold; or even 
in dynamically slowly varying regions of the edge manifold, allowing us to investigate the 
effects of a stable stratification on any coherent structures associated with such states. In 
unstratified plane Couette flow, these coherent structures are manifestations of the self- 
sustaining process (SSP) deduced on physical grounds by Waleffe ( I997| ), or equivalently 
finite Reynolds number solutions of the vortex-wave interaction (VWI) asymptotic equa¬ 
tions initially derived mathematically by Hall & Smith ( |1991 ). The stratified coherent 
states we identify at moderate Reynolds number display an altered form from their un¬ 
stratified counterparts for bulk Richardson numbers Ris = gApH/{prAU‘^) = 0{Re~^), 
and exhibit chaotic motion for larger Rig. We demonstrate that at high Reynolds num¬ 
ber the suppression of vertical motions by stratification strongly disrupts input from 
the waves to the roll velocity structures, thus preventing the waves from reinforcing the 
viscously decaying roll structures adequately, when Rig = 0{Re~^). 


1. Introduction 

Stability theory has been of central importance to fluid dynamics since the pioneering 
work of Reynolds (1883). Linear stability theory is useful for determining whether or 
not for a given situation, all sufficiently small perturbations eventually decay to the 
laminar state. However, it is an inherently nonlinear question to ask what is the domain 
of attraction of a known asymptotically stable solution. The investigation of this second 
question is vital in order to understand why, for example, the laminar solution of plane 
Couette flow (PCF), i.e. the flow between two horizontal plates moving with non-zero 
relative velocity, is linearly stable at every Reynolds number, and yet sustained turbulent 
dynamics have been observed experimentally for Reynolds numbers as low as 325 (see 


Bottin (fe Chate 19981. PCF is therefore a two state system in which both the laminar 
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Turbulence 



0 (energy) 
Laminar state 


Figure 1. Cartoon of state space showing the laminar and turbulent states, the edge manifold, 
edge state and its unstable manifolds, the minimal seed initial condition, and the trajectory of 
the minimal seed through state space. 


state and a chaotic turbulent state are asymptotically attracting, and their respective 
basins of attraction divide the entire state space into two distinct regions, as shown 
schematically in figure These two regions are then separated by an edge manifold, or 
‘edge’. There is substantial evidence to demonstrate that for Reynolds numbers far above 
the critical Reynolds number for transition, below which turbulent dynamics cannot be 
maintained, many shear flows that have an asymptotically attracting laminar state are 
a two state system (see Duguet et al. 20081. Near the critical Reynolds number for 
transition, there is growing evidence (see e.g. Chantry fc Schneider||2014 1 that the edge 
manifold becomes ‘wrapped up’ into the turbulent state, where the precise role of the 
edge manifold is less clear. 

Initial conditions on the edge manifold are attracted to one of a collection of fixed 
points, periodic orbits, relative periodic orbits or chaotic attracting sets residing entirely 
within the edge manifold. The ‘edge state’ is the state in the edge manifold which is 
an attractor for trajectories on the edge. The simplest edge state structure is a single 
saddle in solution space of codimension one whose stable manifold is precisely the edge 
manifold, and whose unstable manifold is directed into the basins of attraction of the 
laminar and turbulent attractors. 

To understand critical conditions for transition to turbulence, it is necessary to identify 
the smallest amplitude perturbation to the laminar state will eventually transition to 
turbulence, which we call the ‘minimal seed’, following Pringle & Kerswell ( |2010[ ). Put 
in the language of edge manifolds, this is the question of what is the ‘closest’ approach 
of the edge manifold to the laminar state, as shown schematically in figure Once 
such a minimal seed is identified, since it lies infinitesimally above the edge manifold, 
and the edge manifold is the stable manifold of the edge state, its trajectory to the 
turbulent attractor will consist of following the edge manifold towards a state in the edge 
manifold before eventually leaving the vicinity of the state along an unstable manifold. 
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towards the turbulent attractor. The edge state and other states in the edge manifold 
therefore determine a key component of the transition mechanism, since the closer an 
initial condition is to the edge manifold, the longer it will spend in the vicinity of such a 
state before being directed towards the turbulent attractor. 


A numerical identification of an approximation to a minimal seed yields valuable sta¬ 
bility information about the laminar state in terms of the minimal energy required for 
transition to turbulence, and the lack of a minimal seed for turbulence corresponds to 
the identihcation of a globally attracting laminar state. Also, the subsequent evolution 
of a sufficiently accurate approximation to a minimal seed should spend a large amount 
of time near a state in the edge manifold and have a long residency time there, before 
moving away from the edge manifold. 


The minimal seed in unstratified PCF found by Rabin et al. (2012) exhibits just such 
a long residency time near a coherent state consisting of nearly streamwise independent 
rolls and streaks, before eventually transitioning to turbulence. This coherent state can 
clearly be interpreted as a manifestation of the ‘self-sustaining process’ (SSP) physically 
described by Waleffe ( 1997| , or equivalently as a finite Reynolds number realisation of 
the asymptotic vortex-wave interaction (VWI) of Hall & Smith (1991) as demonstrated 
numerically by Hall & Sherwin (2010). These solutions to the Navier-Stokes equations 
consist of small amplitude streamwise independent roll structures creating much larger 
amplitude velocity streaks in the flow, which in turn suffer instabilities to produce waves, 
which then reinforce the rolls, hence leading to a self-sustained process. The small ampli¬ 
tude roll structures are readily susceptible to viscous decay, and so the role of the waves 
is to inject energy into the rolls at a rate that exactly offsets this viscous decay. The 
SSP/VWI states are dynamically most sensitive to the feedback from the waves into the 
rolls since this is an inherently nonlinear process, whereas the rest of the cycle depends 
on linear transient growth and linear instability. 


Here, we identify minimal seeds for turbulence in stratified PCF. The addition of a 
stable density stratihcation through hxing the density of the fluid at each of the hori¬ 
zontal plates at different (statically stable) values tends to inhibit vertical motions due 
to them being less energetically favourable. We investigate how the minimal seed and its 
subsequent trajectory in state space, and any identifiable coherent structures on this tra¬ 
jectory, are affected by stratification and its inhibition of vertical motions. In particular, 
we focus on how the energy of the minimal seeds vary with increasing stratihcation, and 
how the coherent states, which in the unstratihed case are SSP/VWI states, vary with 
increasing stratihcation and are disrupted due to the inhibition of vertical motions by 
static stability. 


In section 2 we outline the algorithm used to hnd minimal seeds. We use the same 
nonlinear direct-adjoint-looping (DAL) method used by Rabin et al. (2012), which is the 
nonlinear culmination of optimal perturbation analysis developed hrst for steady linear 


problems, later for time varying linear problems (see Schmid 

2007), and most recently 

for fully nonlinear dynamics (see e.g. 

Pringle & Kerswell 

(2010 

);|Cherubini et al. 

(2010), 

and for reviews 

Luchini & Bottaro ( 

2014) 

Kerswell et al. 

(2014[)). In section 3 we show 


the results of the DAL method applied to stratihed plane Couette how. In section 4 we 
interpret the identihed trajectories of minimal seeds for turbulence in stratihed plane 
Couette how in the language of SSP/VWI states, and demonstrate that the presence of a 
surprisingly weak stable stratihcation can still signihcantly modify the whole interaction 
by disrupting the nonlinear feedback from the waves into the roll structures, principally 
through an inhibition of vertical motions. We draw our conclusions in section 5. 
















































4 


T. S. Eaves and C. P. Caulfield 

2. Direct-Adjoint-Looping (DAL) method 

We consider stably stratified Boussinesq PCF with a linear equation of state in which 
fluid flows between two parallel horizontal plates moving in opposite directions with 
relative speed 2AU, separated by a distance 2H, across which a constant density dif¬ 
ference 2Ap is maintained from a reference density pr 3> Ap. Nondimensionalising with 
respect to At/, H and Ap, and decomposing the total velocity and density fields as 
{utot,Ptot) = [U,p) -I- (u,p), where {U,p) is the laminar state, we obtain 

dfU + [u + U) • {u -\- U) = —Vp — Rispy +Re~^\/‘^u, {u, p){y = ±1) = (0,0) 
dtp -I- (it -I- [/) • V(p -f p) = {RePr)~^V^p, V • it = 0, U = yx, p = —y. (2.1) 


with the three nondimensional parameters: Reynolds number i?e; Prandtl number Pr] 
and bulk Richardon number Ris defined as 


AUH V gApH 

Re = -, Pr = —, Ris = 

V K 


ATT2> ( 2 . 2 ) 

PrAU^ 

where v is the kinematic viscosity and k is the diffusivity of density. For the calculations 


presented here, following Rabin et al. (2012) we choose Re = 1000, set Pr = 1, and vary 
Ris- 

We define a nonlinear optimal perturbation as an initial condition {u,p){t = 0) = 
(uo,po) of given initial total energy (kinetic energy plus potential energy) Bq, defined as 


Eg = (uq ■ uq + RisPo) l‘2-, where {a,b) = ^ fa - bdV, 

y Jv 


(2.3) 


that maximises, over a given time horizon T, a given quantity of interest p,T). 
Here, V is the volume of the computational domain. In order to find minimal seeds 
for turbulence, which are initial conditions whose trajectories eventually transition to 
turbulence, we typically choose an appropriately large large value T = 300, and we 
consider a functional J' that takes heightened values in the turbulent state. With the 
presence of a density field, it is possible to find large amplitude waves that have an 
instantaneously large total energy, but are not a turbulent state. Therefore, rather than 


choosing the total energy density at T as as was done by Rabin et al. (2012) for 


unstratified minimal seeds, we choose the total time averaged dissipation of total energy 
density, the stratified generalisation of the objective functional chosen by |Monokrousos| 
et al. (2011), for an equivalent minimal seed calculation. 


We thus write 


J = {TRe)-^ [Vu : Vu + RisPr-^Vp • Vp] , (2.4) 

where [a, b] = fj' (a, b) dt, and the maximisation of // can be conducted by taking 
variations of the augmented and constrained functional C: 

C = J{u,p,T) - [dtu + N{u) + Vp -I- Rispy - v] - [V • u, q] 

- [dtp + {u + U) • V{p + p) - (PePr)“^V^p, p] -|- {uq - u{0), Vq) 

+ (po - p(0), po) - ((mq, Uq) /2 (Rispl) l‘2. - Eg) c, (2.5) 


where N{u) = (it -f [/) • V(m -1-17). The Lagrange multipliers v, q and p are termed the 
adjoint velocity, pressure and density and together enforce the Boussinesq Navier-Stokes 
equations (2.1) on u and p. The initial conditions and initial energy are enforced by Vg, 
Po and c. Taking variations with respect to u, p, Ug and pg yields the following system of 
‘adjoint’ or ‘dual’ equations that must also be satisfied at all times and points in space 
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Ris 




E 

c (N) 




0 

2.225 

X 

10"' 

3 < 

Ec 

< 

2.250 

X 

10"' 

1 

0 

2.250 

X 

10"' 

® < 

Ec 

< 

2.275 

X 

10"' 

10"® 

2.600 

X 

10"' 

® < 

Ec 

< 

2.625 

X 

10"' 

3 X 10"® 

3.450 

X 

10"' 

® < 

Ec 

< 

3.475 

X 

10"' 

10"® 

6.300 

X 

10"' 

® < 

Ec 

< 

6.400 

X 

10"' 


E, (W) 

8.925 X 10"'^ <Ec< 8.950 x 10”^ 
N/A 
N/A 

1.450 X 10"® <Ec< 1.475 x 10"® 

2.450 X 10"® <Ec< 2.575 x 10"® 


Table 1. Values of the critical energy density Ec, the energy of the minimal seed, for various 
bulk Richardson numbers Rin in the two geometries N and W. The upper bound corresponds 
to the flow evolutions shown in subsequent figures, and in the supplementary materials. The 
lower bound corresponds to an Eq at which a turbulent state cannot be attained. 


by a nonlinear optimal perturbation: 


dtv + {v, u) + Re ^ 
dtp + {U + u) • Vr] 


+ Vq — riV{p + p) — (ReT) = 0 


(RePr) ^W^rj — RiBy 


V • u = 0, 

V - RiB{RePrT)-^V'^p = 0, (2.6) 
v{T) = 0, r 7 (T) = 0, (2.7) 
ry(0) - cpoRiB = 0, (mq, mq) + {R^bpI) - = 0, ( 2 . 8 ) 

where lV^(ui, u) = dj {{Uj + Uj)vi) — Vjdi{Uj + Uj). 

The first step of the DAL method to find a nonlinear optimal perturbation is to choose 
an initial condition guess (tto,po) of energy density Eq. This initial condition is then 
integrated forwards in time to t = T using the direct Boussinesq Navier-Stokes equations 


i ;( 0 ) — cuq = 0 , 


(2.1). These flow fields are stored, and used to integrate backwards in time the adjoint 
fields V and ry from the null ‘end’ conditions (|2.7|) using the adjoint equations (2.6), 


which actually depend on the direct fields {u,p). Once values for the adjoint variables 


are obtained at t = 0, we have compatibility conditions (2.8) relating t)(0) to Uq and 


77 ( 0 ) to po that must be satisfied by a nonlinear optimal perturbation. If not satisfied, 
these compatibility conditions yield gradient information for the objective functional ff 
with respect to changes in the initial condition (uq,po), allowing a refined guess for the 
nonlinear optimal perturbation to be made. This sequence of steps is continued until 
convergence. 

For large T, any initial condition that eventually becomes turbulent will be a turbulent 
seed using this method. The minimal seed, however, has the special property of having 
the lowest critical initial energy density Eq = Ec of all such seeds. We first identify a 
perturbation that causes turbulence, with initial energy density Eq = Ec- This 

initial condition is then used as an initial guess for the DAL method at a smaller initial 
energy density Eq = Ep < E^, with uniformly rescaled energy. The DAL method then 
finds a more efficient route to turbulence at Eq = Ep. This more efficient initial condition 
is then uniformly rescaled to have energy density Eq = E~^ < Ep, and the process is 
repeated. This ‘laddering down’ is continued until Eq = Ec, at which point any further 
reduction in Eq cannot find an initial condition leading to turbulence. 


3. Stratified minimal seeds for turbulence 


We investigate two geometries, namely the ‘narrow’ geometry ‘N’ investigated by|Rabiii 


\et al. (2012) in the unstratified case, 4.357r x 2 x I.OStt, one of the geometries considered 
by Butler & Farrell (1992) for linear, unstratified optimal perturbations, in which the 
unstratified minimal seeds are localised in the streamwise direction but essentially fill the 
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Figure 2. Time variation of energy density E{t) as defined in (3.11 for the minimal seed tra- 

red, 10"" 


jectories in geometry N for Ris =0, black, 10 
pnrple. 


blue, 3 X 10 green and 10 ^ 


spanwise width, and 4.357r x 2 x 2.107r a ‘wide’ geometry ‘W’, which is twice as wide, and 
allows for localisation in the spanwise direction also. We used a modified version of the 
parallelised CFD solver Diablo ( Taylor|2008 ) to solve the forward and adjoint equations, 
which uses Fourier modes in the streamwise x and spanwise z directions, and finite 
differences in the wall-normal y direction, and uses a combined implicit-explicit Runge- 
Kutta-Wray Crank-Nicholson time integration scheme. The resolution for geometry N 
was 128 X 256 x 32 and for geometry W was 128 x 256 x 64. 

Using the laddering down approach described above, we have converged to the minimal 


seed for Ris = 0 (confirming quantitatively the unstratified case of Rabin et al. (2012) 
using a different code, and objective functional J), Ris = 10“^, 10“'^, 3 x 10“'^ and 10“^ 
in geometry N and Ris = 0, 3 x 10“^ and 10“^ in geometry W, using the fixed value 
T = 300 for all bulk Richardson numbers except for the largest, RIb = 10“^, for which we 
use T = 400. The respective values of EfiRis) are shown in table[^ We immediately see 
that Ec is an increasing function of Rib, as expected, since a stable stratification inhibits 
vertical motions, and so a transition process involving vertical motions should be expected 
to require a larger energy input. Interestingly, E,. in geometry W is approximately 40% of 
Ec in geometry N for the same Rib- Since Ec, as defined in (2.3), is an energy density, and 


the volume of geometry W is twice that of geometry N, this suggests that the minimal 
seeds in geometry W are spanwise localised, and that the narrow geometry N actually 
requires higher maximum amplitudes of perturbation in the minimal seed due to the 
enforced spanwise periodicity. 

3.1. Minimal seeds in narrow geometry N 

Figure shows the time evolution from the minimal seed initial conditions for geometry 
N of the total energy density E{t) defined as 


^ “ "b RibP^) = K{t) + P{t), 


(3.1) 


where K{t) is the kinetic energy, P{t) is the potential energy, and angled brackets denote 
volume averaging as defined in ( |2.3[ ). Streamwise velocity u = ±0.6max(u) isosurfaces 
are plotted in figures and at times t = 0 and 25, t = 70 and 150, and t = 210 
and 280 respectively, and videos of the flow evolution are available as supplementary 
materials. 

The unstratified minimal seed, as shown in the leftmost column, consists of a localised 
patch of flow structures aligned against the mean shear, unwrapping via the Orr mech¬ 


anism (see Orr 1907) into an array of streamwise aligned structures with a distinct 


oblique component, shown in the bottom left panel of figure These oblique struc- 























Disruption of SSP/VWI states by a stable stratifieation 


7 



Figure 3. Isosurfaces of streamwise perturbation velocity u = ±0.6max(u) at t = 0 (top 
row) and t = 25 (bottom row) for the minimal seed trajectories in geometry N for Ris = 0, 
10■^ 3 X 10 and 10 from left to right. Videos of the flow evolution are available as 
supplementary materials. 


tures then transfer energy into streamwise independent streaks through the oblique wave 
mechanism in which structures with wavenumbers {k^, kz) = (0, ±o) interact nonlinearly 
to move energy into the wavenumeber {kx,kz) = (0,0). These streaks are then able to 
‘self-sustain’ for a long period of the flow’s evolution, as shown in the leftmost column of 
figure]^ by utilising the lift-up mechanism (described by Landahl 1980) and their own 
instability to offset viscous decay, (see Waleffe 1997), before eventually being of large 
enough amplitude to transition to a high energy oblique structure, as shown in the left¬ 
most column of figure which is visited only transiently, leading to a break down to 
small scale turbulence due to an instability reminiscent of the Kelvin-Helmholtz insta¬ 
bility. The energetics of these sequential growth mechanisms in unstratified minimal seed 


trajectories were discussed by 

Pringle et al. 

(2012 

) and 

Duguet et al. 

(20K 

il 

This transition mechanism has been interpreted by 

Rabin et al. ( 

2012 ) 

a 


concentration of energy in a small region of the flow, allowing a minimisation of the 
total energy input whilst maximising the local energy input. From this point, the most 
efficient way to transition to turbulence is by exploiting the lift-up mechanism on a high 
energy flow structure and thus increasing the total energy of the flow to a point beyond 
which it is sufficiently unstable to become turbulent. As shown in figure there is a 
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Figure 4. Isosurfaces of streamwise perturbation velocity u = ±0.6max(u) at t = 70 (top 
row) and t = 150 (bottom row) for the minimal seed trajectories in geometry N for Ris = 0, 
10~^, 10“^, 3 X 10“^, and 10“^, from left to right. Videos of the flow evolntion are available as 
snpplementary materials. 


sustained period of near-constant energy denstity E{t), and during this time the flow is 
only very slowly changing, indicating that the flow trajectory is near a coherent state. 
The isosurfaces plotted in the leftmost columns of hgure show that this state is very 
reminiscent of an SSP/VWI state. As Ris is increased to Ris = 10“'^ and 10“^, the 
minimal seed, its trajectory, and the coherent state to which it approaches, remain very 
similar to the unstratified case, as can be seen from the flow isosurfaces shown in the 
second and third columns of figures and although there are the beginnings of an 
oblique structure appearing in the coherent states to which these flows evolve. 

As Rib is increased further, the behaviour changes qualitatively. The minimal seed 
for Rib = 3 x 10“^, (as shown in the fourth column, top row of figure still consists 
of a localised patch of flow structures aligned against the mean shear, which unwrap 
via the Orr mechanism into the same array of streamwise aligned structures with a 
distinct oblique component, followed by the oblique wave mechanism. However, the newly 
created streamwise independent streaks are visited only transiently, no longer able to be 
sustained, and the flow evolves into a new, fully three dimensional coherent state which is 
shown in the fourth column of figure]^ This new stratified coherent state is different from 
the unstratified SSP/VWI state visited by the unstratified minimal seed trajectory, with 









Disruption of SSP/VWI states by a stable stratifieation 


9 



Figure 5. Isosurfaces of streamwise perturbation velocity u = ±0.6max(u) at t = 210 (top 
row) and t = 280 (bottom row) for the minimal seed trajectories in geometry N for Ris = 0, 
lO"'^, 10■^ 3 X 10 and 10 from left to right. Videos of the flow evolution are available as 
supplementary materials. 


the stratified coherent state being clearly more three dimensional and very reminiscent 


20141. 


of the stratified edge states recently reported by Olvera & Kerswell 

Figure shows that in the cases Ris = 10“'^, 10“^ and 3 x 10 transition to tur¬ 
bulence is apparently faster than the unstratified case, with the same accuracy in the 
estimated value of i?c, and the same time target time T = 300. However, the time taken 
to transition to turbulence is a function of how close our numerical estimates for the 
minimal seed initial conditions are to the edge manifold. In the limit of successively bet¬ 
ter approximations to the ‘ideal’ minimal seed, which lies exactly on the edge manifold, 
Uransition —t OO- Therefore, since the value of Ec is only bracketed to a certain accuracy, 
the time to transition for the minimal seeds presented here is most likely to be affected 
by the actual value of relative to the energy bracket for Ec found here. 

The minimal seed for Rig = 10“^ is once again of a qualitatively different character. 
Although it still consists of a localised patch of flow structures that unwrap via the Orr 
mechanism, there appears to be no quasi-steady flow structure into which it evolves. 
The flow is now chaotic with a weak oscillation, with no single structure dominating the 
flow evolution. Indeed, due to the chaotic nature of the dynamics on the edge manifold 
that this trajectory follows, the DAL method struggled to identify turbulence solutions 
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without extending the optimisation time interval to T = 400. Even after extending the 
optimisation time interval, the DAL method required up to ten times as many iterations 
to identify initial conditions that transition to turbulence for the flow with Ris = 10“^ 
compared to the number of iterations required for flows with smaller bulk Richardson 
numbers, apparently because of the chaotic nature of the edge manifold. As noted in the 
hgure captions, videos of the evolution of the minimal seeds for the different values of 
Rib are available as supplementary material. 


3.2. ‘Wide’ geometry W 

As already noted, the minimal seeds in the narrow geometry N are streamwise localized, 
but fill much of the spanwise extent of the computational domain. To investigate to 
some extent the sensitivity of the identified minimal seeds to the flow geometry, we also 
calculate minimal seeds in geometry W, which is twice as wide in the spanwise direction. 
Figure shows the time evolution from the minimal seed initial conditions of the total 
energy density E{t) as defined in (3.11 for the wide geometry W. Streamwise velocity 
u = ±0.6 max(u) isosurfaces are plotted in hgures and at times t = 0 and 25, 
t = 70 and 150, and t = 210 and 280 respectively, and videos of the flow evolution are 
available as supplementary materials. 

The RiB = 0 unstratified minimal seed trajectory in this wider geometry shares the 
same characteristic evolution as the minimal seed trajectory in the narrower geometry, 
but with the addition of spanwise localisation, that is a streamwise and spanwise localised 
patch of flow structures aligned against the mean shear which unwrap via the Orr mech¬ 
anism into a streamwise aligned structure with a distinct oblique component, as shown in 
the leftmost column of figure before utilising the oblique wave mechanism to produce 
a spanwise isolated pair of streaks, reminiscent of the non-localised structure seen in ge¬ 
ometry N. As is apparent in the tabulated values of E,, listed in table these structures 
are slightly less energetic than the minimal seeds identihed in the narrow geometry N, 
in that, as already noted, the critical values of the energy density Ec (i.e. the energy 
divided by the volume of the computational domain) in geometry W are approximately 
40% of the equivalent values determined in geometry N. 

These streaks survive in the flow for an extended period of time, as shown in the left¬ 
most column of figure and are another realisation of an SSP/VWI coherent structure. 
These streaks eventually transition to a high energy spanwise isolated oblique structure, 
as shown in the leftmost column of figure which is visited only transiently, before 
breaking down to small scale turbulence. This sequence of events, as well as the span- 
wise localisation, are both consistent with those reported by Monokrousos et al. (2011) 
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13.66 
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13.66 




6.62 


13.66 


Figure 7. Isosurfaces of streamwise perturbation velocity u = ±0.6 max(u) at t = 0 (top row) 
and t = 25 (bottom row) for the minimal seed trajectories in geometry W for Rib = 0, 3 x 10“^, 
and 10~^, from left to right. Videos of the flow evolution are available as supplementary materials 


and verified using a different objective functional by Rabin et al. (2012) in the domain 
Itt X 2 X 27 r at the larger Reynolds number Re = 1500. 

The minimal seed trajectory for = 3 x 10“^ (shown in the middle column of figures 
[^[^an again shares the same characteristic evolution as the equivalent minimal seed 
trajectory in the narrow flow geometry N. The initial condition consists of the same 
spanwise and streamwise localised patch of flow structures aligned against the mean 
shear as in the unstratified flow, which again unwrap via the Orr mechanism into a 
streamwise aligned structure with a distinct oblique component, as shown in the lower 
row, middle column of figure The oblique wave mechanism transfers energy into a pair 
of spanwise isolated streaks that are visited only transiently, and the flow evolves onto a 
long lived spanwise localised three dimensional coherent state which has the same oblique 
characteristics as the one found in the narrow geometry N. This structure is eventually 
no longer able to be maintained, and transition to turbulence occurs. 

Furthermore, the minimal seed trajectory for RIb = 10“^ in geometry W is also a 
spanwise localised version of the equivalent trajectory in the narrower geometry N, as 
can be seen by comparison of the rightmost columns of figures and and figures [7j|^ 
an m The initial condition in geometry W consists of a streamwise and spanwise localised 
patch of flow structures that unwrap via the Orr mechanism into a flow which quickly 
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Figure 8. Isosurfaces of streamwise perturbation velocity u = ±0.6max(u) at t = 70 (top 
row) and t — 150 (bottom row) for the minimal seed trajectories in geometry W for Ris = 0, 
3 X 10“®, and 10“^, from left to right. Videos of the flow evolution are available as supplementary 
materials. 


becomes chaotic. The flow continues in much the same way as its small domain version, 
until eventually transitioning to turbulence. Once again, it appeared computationally 
more difficult to identify the minimal seed for the flow with Ris = 10“^ than for the flows 
with smaller bulk Richardson numbers, requiring two or three times as many iterations, 
although this convergence was still faster than for the equivalent flow with the same 
Rib = 10“^ in the narrower geometry N. 


4. Rolls, streaks and waves 

To analyse the effect that stratification has on the SSP/VWI states, we decompose 
the full perturbation velocity and density fields into roll, streak and wave components. 
First, we define U = {u)^ and 0 = {p)^, where {a)^ = a da;, and decompose 

u = U + u = Ur + Us + u and p = 0 + p, so that 

(^U, = (hir pUs a, p) = (0,V, W , o') roll T (Id ,0,0, Q') streak “t” p)wave-> (4-1) 


where the r subscript denotes the streamwise-independent wall-normal and spanwise 
roll velocity and the s subscript denotes the streamwise-independent streamwise streak 
velocity. As originally argued independently bylHall & Smithl (|1991|) andlWaleffel (|1997|), 
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Figure 9. Isosurfaces of streamwise perturbation velocity u = ±0.6max(u) at t = 210 (top 
row) and t — 280 (bottom row) for the minimal seed trajectories in geometry W for Ris = 0, 
3 X 10“®, and 10“^, from left to right. Videos of the flow evolution are available as supplementary 
materials. 


at high Reynolds number, if there are rolls in the flow of typical amplitude 0(e) with 
e< 1, then their decay rate due to viscosity is 0{Re ^). During the 0{Re) time in which 
they survive in the flow, they can advect streamwise velocity through the 0(1) shear of 
PCF a distance O(eRe) and so produce 0{eRe) streaks in the flow. If the amplitude of 
these streaks is sufficiently large, they can undergo an instability which creates a wave 
field. The nonlinear self interaction of this wave field then puts energy back into the rolls, 
and provided that this input of energy is sufficient to balance the viscous decay of the 
rolls, a ‘self-sustaining process’ associated with this ‘vortex-wave interaction’ is possible. 
Supposing that the streaks need to be 0(1) to become unstable requires e = Re~^ 
and so a quadratic self-interaction of the wave field of 0{Re~^) is sufficient. Hall & 


Sherwin (2010) showed that there is a numerically exact solution in unstratified PCF for 
asymptotically large Re such that the rolls have V,W = 0{Re~^) and the streaks have 
U = 0(1) throughout the domain. The waves inject energy into the rolls primarily in an 
0{Re~^^^) critical layer where the background PCF plus the streak flow has the same 
velocity as the wave velocity, within which (u, v, w) = 0(i?e“®^®, Re~'^^^, Re~^^^), giving 
an integrated Reynolds stress contribution over the critical layer of 0{Re~^). The waves 
essentially provide an 0{Re~^) jump across the critical layer in the component of the 
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Figu re 10. Time dependence of normalised streak kinetic energy density KsjK as defined in 
( |4-3|) (p lotted with a blue line) and normalised wave kinetic energy density K^/K as defined 
in ( |4.5 1 (plotted with a red line) for the minimal seed trajectories in geometry N for Rib ~ 0 
(upper left), 10“'* (upper middle), 10“® (upper right), 3 x 10“® (lower left), 10“^ (lower right). 
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Figu re 11. Time dependence of normalised streak kinetic energy density Kg/K as dehned in 
( |4.3|) (p lotted with a blue line) and normalised wave kinetic energy density K^/K as defined 
in ( |4.5 1 (plotted with a red line) for the minimal seed trajectories in geometry W for Rib = 0 
(left), 3 X 10“^ (middle) and 10“^ (right) 


roll velocity normal to the critical layer through their Reynolds stresses. Such scalings 
have been used successfully to provide initial guesses for a search for numerically exact 
solutions of the Navier-Stokes equations (see WaleffepOOl Wedin &: Kersweli||2004 ). 

To see how this scaling is affected by the addition of a stable stratification, we examine 
the energetics for the total energy density of the rolls and of the streaks. We obtain 


^Kr{t) = \^(\Ur\’^) ={-RibVQ-U 
dt 2 dt ' 


ii ■ Vii)^ 

-Re-^VUr : , 

-Re-^VUg-.VUg)^^^, 

\ - RiB<d {u • Vp)^ 

-RiB{RePr)~^Ve-Ve) 


Aif (i) = lA|2\ 
dt ^ ^ 2dt '' ' ' 




(4.2) 


(4.3) 


(4.4) 


where ^ a dzdy, Kr{t) is the roll kinetic energy density, Kgit) is the 

streak kinetic energy density, and Ps{t) is the streak potential energy density. Note that 
Kfit) + Ks{t) + Kyj{t) = K{t), where K^{t) is the wave kinetic energy density 


K^{t) = -{\vi\^) 


(4.5) 


Taking e = Re th e correct scaling for the viscous decay time is obtained from ( |4.2[ ) . 
We also see from (4.3) that when V = 0(i?e“*) is sustained over a period of 0{Re), we 
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obtain streaks U = 0(1). We can also obtain the scaling of Hall & Sherwin (2010) for 
the wave velocities, and their Reynolds stress contribution, under the assumption that 
they act over an 0{Re~^^^) critical layer, and balance the viscous dissipation there. 

For stratified PCF with Ris 0 there is a new buoyancy flux term —RisVQ entering 
the energetics of the rolls. Comparing the first term in the right-hand side of (4.3) to 
that of (4.4), which are production terms for the streak kinetic energy density and streak 
potential energy density respectively, we see that 0 = 0(1) everywhere for a passive 
scalar field placed in a self-sustained process. Thus, (4.2) shows that a buoyancy flux 
associated with the streak flow effective across the whole domain provides a contribution 
to the kinetic energy density of the rolls of 0{RiBRe~^), and is able to disrupt fully 
the flux of wave Reynolds stress input of 0{Re~^) over the 0{Re~^^^) critical layer when 
Rib = 0{Re~^). This simple minded high Reynolds number scaling argument, which is 
domain size independent, is not inconsistent with the moderate Reynolds number minimal 
seed calculations presented above for which the unstratified coherent state visited by the 
unstratified minimal seed trajectory is no longer a viable solution for sufficiently large 
Ris ~ 3 X 10“^, and the new coherent states have a modified roll structure that is no 
longer streamwise-independent. 

Figures 10 and 11 show the time evolution of the streak kinetic energy density Ks{t), as 
defined in (4.3) and plotted with a blue line, and the wave kinetic energy density K^[t), 
as defined in (4.5) and plotted with a red line, normalised by the total kinetic energy 
K{t) for each of the minimal seed trajectories in geometries N and W respectively. It is 
clear that for RiB =0, 10“'* and 10“^, there is a balance (shown by the approximate 
plateaux in the streak and wave energy components) between these energies (and also 
the residual roll kinetic density density Kr{t), which is not shown as it is, unsurprisingly, 
appreciably smaller in magnitude) for a large period of the flow evolution. For RiB = 
3 X 10“^, this balance has been significantly disrupted, and is not maintained purely 
by the velocity fields. There also is a small amplitude oscillation, and for RiB = 10“^ 
the above-mentioned weakly oscillatory nature of the flow is clear. The unstratified self- 
sustaining process is completely disrupted for Rfs > 3 x 10“^ as expected. 


5. Discussion 

Using the direct-adjoint looping (DAL) method, we have computed minimal seeds for 
turbulence (the initial conditions of smallest possible initial perturbation energy density 
Ec{RiB) that transition to turbulence) in stratified PCF for a range of bulk Richardson 
numbers RiB in two geometries: a narrow geometry labelled N, which apparently allows 
only streamwise localisation of the minimal seed initial condition; and a twice as wide 
geometry W which also allows spanwise localisation. In the unstratified case we converge 


to the same minimal seed found by Rabin et al. (2012) in geometry N, although here we 


use the time averaged dissipation rate objective functional in the DAL method instead 
of the total energy density at the target time, and a minimal seed very similar to that of 


Monokrousos et al. (2011) in geometry W. 


Since a stable stratification inhibits vertical motions, we see an increase of Ec with 
Rib, as expected. The minimal seeds follow trajectories in state space close to the edge 
manifold towards a state in the edge manifold. These trajectories are found to spend a 
large amount of time in the vicinity of such a state. For unstratified flows, these coherent 
states are a realisation of an SSP/VWI state, and for sufficiently small RIb the coherent 
states are largely unchanged. For RIb = 0{Re~^) there is a slow draining of energy into 
the density field, creating a highly modified, yet still stationary coherent state, while for 
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larger Ri b the coherent states feature inherently three-dimensional chaotic motion with 
weak oscillations. 

Examining the flow in terms of roll, streak and wave components as defined in (4.1) 
demonstrates that the density field is expected, at asymptotically higher Reynolds num¬ 
ber, to have a significant disrupting effect on the SSP/VWI process, and hence the 
coherent states, when Rib = 0{Re~^) through a removal of energy from vertical motion 
in the rolls, which is not entirely inconsistent with the minimal seed trajectories found 
at moderate Reynolds number. The effects of stratification on a wide class of exact co¬ 
herent states in shear flows have yet to be studied in detail. We have demonstrated here 
that stratification disrupts the well-established high Reynolds number SSP/VWI states 
in PCF for the small value Rib = 0{Re~^) by affecting the energy input into the roll 
structures, the most delicate part of the interaction process, through an inhibition of 
vertical motions. For larger bulk Richardson numbers, we also observe chaotic solutions. 
We believe therefore that a careful re-examination of the SSP/VWI ansatz must be made 
for the case of stratified shear flows even with a very weak stratification when i?e ^ 1, 
an important class of flows common in both environmental and industrial contexts. 
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